Radiative energy loss in a finite dynamical QCD medium 
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The radiative energy loss of a quark jet traversing a finite size QCD medium with dynamical 
constituents is calculated to first order in opacity. Although finite size corrections reduce the energy 
loss relative to an infinite dynamical QCD medium, under realistic conditions it remains significantly 
larger than in a static medium. Quantitative predictions of jet suppression in relativistic heavy ion 
collisions must therefore account for the dynamics of the medium's constituents. Finite size effects 
are shown to induce a non-linear path length dependence of the energy loss. Our results suggest a 
simple general mapping between energy loss expressions for static and dynamical QCD media. 
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1. Studying the suppression of high transverse momen- 
tum hadrons is a powerful tool to map out the density of 
a QCD plasma created in ultrarelativistic heavy ion colli- 
sions P, 0]. Since this suppression (called jet quenching) 
results from energy loss of fast partons moving through 
the plasma quantitative jet quenching predic- 
tions require reliable energy loss calculations. 

In the majority of currently available studies the 
medium-induced radiative energy loss is computed by 
assuming that the QCD medium consists of ran- 
domly distributed static scattering centers ("static QCD 
medium"). We recently calculated Q, at leading order 
in opacity, the heavy quark radiative energy loss in an in- 
finite QCD medium consisting of dynamical constituents 
and found that the energy loss increases by almost a fac- 
tor 2 relative to an equally dense static medium. How- 
ever, this calculation was performed in the Bcthc-Hcitler 
limit which is well known 0] to overpredict radiative en- 
ergy loss since it does not include coherence and finite 
size effects. As the medium created in heavy ion col- 
lisions has finite size, it is essential to explore how the 
qualitative conclusions obtained in Q change once such 
effects are included. We find that finite size and coher- 
ence effects decrease the radiative energy loss per unit 
path length more strongly in a dynamical than in a static 
medium, reducing the energy loss ratio between equally 
dense dynamical and static media. Still, the ratio re- 
mains significantly larger than unity even if the medium 
is finite, showing that for quantitative predictions of ra- 
diative energy loss it is important to account for the dy- 
namic nature of the QCD medium's constituents. 

2. We briefly outline the computation of the medium 
induced radiative energy loss for a heavy quark to first 
order in opacity. We consider a QCD medium of size L 
and assume that the heavy quark is produced at time 
Xq = at the left edge of the medium, traveling right. 
Collisions with partons in the medium induce the radi- 
ation of gluons, causing the quark to lose energy. The 
radiative energy loss rate can be expanded in the number 
of scattering events suffered by the heavy quark. This is 



equivalent to an expansion in powers of the opacity. For 
a finite medium, the opacity is given by the product of 
the medium density with the transport cross section, in- 
tegrated along the path of the heavy quark. The leading 
(first order) contribution corresponds to one collisional 
interaction with the medium, accompanied by the emis- 
sion of a single gluon. This is the process we compute. 

To introduce the finite size of the medium we start from 
the approach described in Q and follow the procedure 
used in (loj . The medium extends for a length L from the 
production point of the energetic heavy quark, and the 
collisional interaction inducing the radiation of a gluon 
occurs after a distance I < L inside the medium. 

As in 0, we describe the medium by a thermalized 
quark-gluon plasma at temperature T and zero baryon 
density, with n / effective massless quark flavors in equi- 
librium with the gluons. Three typical Fcynman dia- 
grams contributing to the radiative quark energy loss 
at first order in opacity are shown in Fig. [TJ The di- 
agrams are evaluated in finite temperature field theory 



111 . Il2| , using Hard Thermal Loop (HTL) resumed prop- 
agators [IH for all gluons. A full account of the calcu- 
lation will be presented elsewhere [l3| ■ A flavor of what 
it involves is given by Fig. [TJ The elliptic blob repre- 
sents a source J which at time a;o produces an energetic 
quark with momentum p'. In diagram [JJA_ the produced 
quark is on-shell. It first radiates a gluon with momen- 
tum k = (uj,k z ,k) and then exchanges a virtual gluon of 
momentum q = (qo,q z , <z) with a parton in the medium, 
finally emerging (at the dashed line denoting an on-shell 
cut through the amplitude represented by the diagram) 
with (measured) momentum p = (E,p z ,p) [14j . Since 
the energetic quark produced by the source J can be off- 
shell, we also have contributions such as those in Figs.[TjB 
and[Tp. Amplitudes [TJB and[Tp interfere with amplitude 
[TJ\, leading to the appearance of LPM-like effects once 
all relevant contributions are summed. The present cal- 
culation differs from that in Q by the use of HTL gluon 
propagators to describe the interaction of the quark with 
the medium, and from that in Rcf. by allowing the 
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FIG. 1: Three typical 
Feynman diagrams con- 
tributing to the quark ra- 
diative energy loss in a fi- 
nite size dynamical QCD 
medium at first order in 
opacity. See text for 
discussion. 



jet to be on- or off-shell and restricting the vertices cor- 
responding to gluon exchange to be located inside the 
medium, i.e. at I < L. We use the same kinematic 
approximations as in @, @] ! accordingly, the gluon prop- 
agators for exchanged gluons in Fig. [T] contribute only 
for space-like momenta (qo < |q|) and those for radiated 
gluons only for time-like momenta (u > |k|) @, We 
also assume that J changes slowly, i.e. J{p') ~ J{p) [TBI ] . 

Explicit calculation of all 24 diagrams contributing 
to first order in the opacity yields the follow- 

ing expression for the fractional radiative energy loss 
(/j = gT \jN c /'d + Nf/6 is the Debye screening mass and 
parametrizes the density of the medium): 
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stant, and Cr=|. Further, \ = M 2 x 2 + m 2 where x is 
the longitudinal momentum fraction of the heavy quark 
carried away by the emitted gluon and m g = is the 
effective mass for gluons with hard momenta k > T. 

Similar to the infinite medium studied in Q, each in- 
dividual diagram contributing to the energy loss in a fi- 
nite dynamical medium diverges logarithmically in the 
limit of zero transverse momentum exchange, q—> fl3j |. 
In a dynamical QCD medium both transverse and lon- 
gitudinal gluon exchange contribute to radiative energy 
loss |16[; while Debye screening renders the longitudinal 
gluon exchange infrared finite, transverse gluon exchange 
causes a logarithmic singularity due to the absence of 
magnetic screening (l2| . Remarkably, this singularity is 
found to cancel in the sum over all diagrams [13| . natu- 
rally regulating the energy loss rate. 

3. We can compare the radiative energy loss rate in 
a dynamical medium ([T]) to the analogous result for a 
static medium. One can rewrite the DGLV expression Q 
for the first order radiative energy loss in a static QGP, 
AE s t a t/E, in the same form as Eq. (fl]), except for two 
simple substitutions: (1) Ad yn is replaced by the "static 
mean free path" A sta t , defined by 0, LL7| 
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where c(nf) = 6 1, ^° 2 pj^-7g is a slowly increasing func- 
tion of nf that varies between c(0) m 0.73 and c(oo) w 
1.09. For n f = 2.5 (see below) c(2.5) w 0.84. (2) The ef- 
fective cross section under the integral ([1]) for the energy 
loss rate is replaced as 
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Taken together, these differences will be seen to cause a 
significant increase of the heavy quark energy loss rate in 
dynamical compared to static QCD media. 

These two simple replacements are identical to those 
found in the Bethe-Heitler limit Q. The simplicity of 
this substitution rule is surprising, given the complexity 
of the calculations and their different structure for static 
Q and dynamical (l3j media. (Remember the infrared 
divergences in the dynamical case which cancel only af- 
ter summing all 24 diagrams but don't arise at all in the 
static case.) The integrands in Eq. (fl]) and its static ana- 
log are significantly different from the corresponding ones 
in the Bethe-Heitler limit giving rise (as we will see) 
to a different energy dependence of the dynamic/static 
energy loss ratio. Nonetheless the same simple substitu- 
tion rule is found to apply, suggesting a possibly general 
mapping between static and dynamic QCD media. 

The study presented here considers a finite, optically 
thin dynamical QCD medium (QGP), extending the 
DGLV approach [8j to include parton recoil. In this sense 
it is c omp lementary to the work by Arnold, Moore and 
Yaffc [18I | who study energy loss in an infinite, optically 
thick QGP. We note that the AMY approach [l|] yields 
the same form ([3]) for the effective cross section in a dy- 
namical QCD medium as found here (see also [3]), sup- 
porting our conjecture above. 

4- We now highlight finite size effects in a dynami- 
cal QCD medium, to first order in opacity, with a few 
numerical results for radiative energy loss. In Fig. [2] 
we show the ratio of the radiative energy loss rates in 
equally dense dynamical and static QCD media as a func- 
tion of the initial energy of the fast quark, under RHIC 
and LHC conditions. In both cases a medium of length 
L = 5 fm, a constant value of a s =0.3, and a chemically 
equilibrated QGP with nf = 2.5 effective light quark fla- 
vors is assumed. The light quark mass is assumed to be 
dominated by the thermal mass, M q = /j,/^/6, where /i is 
the Debye screening mass. For the charm and bottom 
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FIG. 2: Ratio of the radiative 
energy loss in finite dynami- 
cal and static QCD media of 
length L = 5fm for light, 
charm and bottom quarks 
(left, center, and right pan- 
els, respectively), as a func- 
tion of initial quark momen- 
tum p. Top row: RHIC 
conditions (average medium 
temperature T = 225 MeV). 
Bottom row: LHC conditions 
(T = 400 MeV). The dashed 
curves show the corresponding 
energy loss ratio in an infinite 
QCD medium for comparison. 



masses we use M c = 1.2 GeV and Mb = 4.75 GeV, respec- 
tively. For Au+Au collisions at top RHIC energies we as- 
sume an average medium temperature of T = 225 MeV, 
for Pb+Pb at the LHC we take T = 400 MeV. 

In all cases, the energy loss is seen to be significantly 
larger in the dynamical than in the static medium. A 
common factor to all situations (large and small jet quark 
masses, hotter and cooler media, finite and infinite me- 
dia) is the O(20%) increase of the energy loss in dy- 
namical media arising from the shorter mean free path 
Adyn ~ 0.84A sta t- The additional increase arising from 
the change ([3]) in cross section is larger; in the energy 
range shown in Fig. [5] it ranges from about 25% to over 
100%, depending on medium temperature and the mass 
and energy of the fast quark. The reduction of the en- 
ergy loss ratio ^^ dyn by finite size corrections is seen to 
be larger for lighter quarks and larger jet energies. The 
smallest finite size corrections and, in the end, the biggest 
dynamical increase are seen for bottom quarks at RHIC. 

Furthermore, in Q we found that dynamical medium 
effects are largest for light quarks and decrease with quark 
mass. Here we see the opposite tendency: after finite size 



correction the energy loss ratio 



AE r , 



becomes smallest 



for light quarks, increasing with quark mass. Such behav- 
ior is important, since it may contribute toward under- 
standing the observed large suppression of non-photonic 
electrons in central Au+Au collisions at RHIC |21| . 

Figure [3] shows that the strength of the finite size cor- 
rections correlates with the dependence of the fractional 
energy loss on the thickness L of the medium. In the 
Bethe-Hcitlcr limit studied for infinite media in Ref. Q, 
quarks of all masses and energies lose energy at fixed 
rate resulting in a linear dependence of the radia- 



tive energy loss on the length L traveled by the quark. 
In contrast, Fig. [3] shows a non- linear L-dependence that 
becomes perfectly quadratic (correspondin g to the deep 
Landau-Pomeranchuk-Migdal (LPM) limit [20( ) for large 
jet energy (see Eq. (|4]) below). The weakest deviations 
from the linear Bethe-Heitler L-dependence are seen for 
low-energy bottom quarks at RHIC, where Fig. [2] (top 
row, right panel) also shows the smallest finite size cor- 
rection. The i-depcndcncc is closest to quadratic for 
light quarks and for very energetic charm and bottom 
quarks at the LHC where also the finite size effects are 
largest. This shows that the finite size corrections imple- 
mented in the present calculation simulate the destruc- 
tive effects of LPM" interference in an infinite medium 41. 



This behavior is expected 15j since the nuclear medium 
has finite dimensions that may be small compared to the 
jet radiation coherence length, especially in the case of 
light partons or high jet energies. Due to this, in finite 
size media the basic formation time physics developed by 
LPM [2(| leads to strong destructive interference effects 
on the quark quenching, as observed in Fig [3] 

5. We finally point out that, contrary to the first order 
study in the Bethe-Heitler limit Q where the energy loss 
ratio Ag d t Y ° saturates for sufficiently large quark energies, 
this ratio keeps decreasing with increasing quark energy 
once finite size corrections are accounted for. In fact, 
one finds analytically [IH that for asymptotically large 
jet energies Eq. |T]) reduces to 
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FIG. 3: First order fractional radiative energy loss as a function of medium thickness L for initial jet energies E — 10, 25, and 
500 GeV (left, center, and right panels, respectively). The two left panels correspond to RHIC conditions, the right panel to 
LHC conditions. Solid, dashed and dot-dashed lines describe light, charm, and bottom quark energy loss, respectively. 



and that the energy loss ratio approaches 

a zp \ ln±§£ \ 

J™ AF ' = p lm 1 1 iET -i = 1 • ( 5 ) 

E^oo A£/ Sta t B-oo A dyll 111 ^ 1 Adyn 

The static approximation thus becomes valid for asymp- 
totically large jet energies. 

6. In summary, we have presented a calculation to 
first order in opacity of the radiative energy loss of a fast 
quark traveling through a finite dynamical QCD medium. 
Finite size effects are found to be most important in the 
ultrarelativistic limit and they effectively reproduce the 
effects of destructive Landau-Pomeranchuk-Migdal inter- 
ference. The calculation suggests the possibility of a 
general mapping between the energy loss expressions for 
static and dynamical media, which we conjecture to carry 
over to higher-order calculations. It also shows that the 
approximation of the medium by a random distribution 
of static scatterers becomes valid in the limit of asymp- 
totically large jet energies once finite size and LPM in- 
terference effects are taken into account. For realistic jet 
energies and medium temperatures reachable at RHIC 
and LHC, however, parton recoil in the medium must be 
accounted for and leads to a large (40-70%) increase of ra- 
diative energy loss when compared with an equally dense 
static medium. This effect is largest for bottom quarks 
at RHIC which may be important for understanding the 
observed large suppression of non- pho tonic electrons in 
central Au+Au collisions at RHIC [2lT |. 
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